Longitudinal and transverse forces on a vortex in 
superfluid 4 He 



H M Cataldo, M A Desposito, and D M Jezek 

Departamento de Ffsica, Facultad de Ciencias Exactas y Naturales, 

Universidad de Buenos Aires, RA-1428 Buenos Aires, Argentina 

and Consejo Nacional de Investigaciones Cientfficas y Tecnicas, Argentina 

Abstract. 

We examine the transverse and longitudinal components of the drag force upon a 
straight vortex line due to the scattering of liquid 4 He excitations. For this purpose, 
we consider a recently proposed Hamiltonian that describes the dissipative motion 
of a vortex, giving an explicit expression for the vortex-quasiparticle interaction. 
The involved dissipative coefficients are obtained in terms of the reservoir correlation 
function. Most of our explicit calculations are concerned to the range of temperatures 
below 0.4 K, at which the reservoir is composed by phonon quasiparticle excitations. 
We also discuss some important implications in the determination of possible scattering 
processes leading to dissipation, according to the values of vortex mass found in the 
literature. 
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Although it is well established that vortex dynamics plays an important role in 
the behavior of superfluids, there are many controversial results in the literature about 
its dissipative motion. At zero temperature the motion of the vortex is provoked by 
the transversal Magnus force, but at finite temperatures the scattering of collective 
excitations by the vortex provides a dissipation mechanism that damps its cyclotron 
motion. It is widely accepted that for low enough temperatures, and only taking into 
account phonon-vortex scattering, the longitudinal component of this dissipative force 
behaves as fifth power in temperature fl|]. Nevertheless, the form of the transverse 
component is still controversial || since very different results were reported so far and 
it does not appear that any of them can be regarded as definitive. 

In this context, we have recently presented || a model for the dissipation of a 
straight vortex line in superfluid 4 He in which the vortex is regarded as a macroscopic 
quantum particle whose irreversible dynamics can be described in the frame of 
Generalized Master Equations. This procedure enables us to cast the effect of the 
coupling between vortex and heat bath as a drag force with one reactive and one 
dissipative component, in agreement with phenomenological theories. In this paper 
we shall investigate the components of the drag force considering that the reservoir is 
composed by the quasiparticle (qp) excitations of the superfluid. The dissipative vortex 
dynamics arises then from the scattering processes ocurring between the vortex and the 
qp's. Most of our explicit calculations will be concerned to the range of temperatures 
below 0.4 K, at which the qp excitations constitute a phonon reservoir Q]. 

Let us begin by performing a description of our model. Choosing a coordinate 
system fixed to the superfluid, the Hamiltonian for a straight vortex line parallel to the 
z-axis, may be written as [[3], [4J 

H v = [P - <Z„A(r)] 2 , (1) 

where r and p denote respectively the vortex position and momentum operators, 

A(r) = ^(y,-x) (2) 

is the vector potential for the Magnus force, m v is the inertial mass of the vortex, p s the 
number density of the superfluid, h Planck's constant, I the system length along the z— 
axis and q v = ±1 the sign of the vorticity according to the right handed convention. 

In what follows we shall consider the excitations to be at rest, so that the velocity 
of the normal fluid vanishes. In a previous work || we have considered an interaction 
Hamiltonian of the form 

H mt = -B v (3) 

where B and v are vectors that depend on the reservoir and vortex operators 
respectively. We have proven that the only linear combination of the vortex observables 
that leads to a dynamics in accordance with phenomenological descriptions is, 
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being Vt the "cyclotron frequency" , 

n=*W. (5) 

m v 

In order to obtain an equation of motion for the mean value of the complex vortex 
position operator R = x + iy, in previous works [[3], ||] we have derived by means of a 
standard reduction-projection procedure, a generalized master equation for the density 
operator of the vortex. We have employed a usual weak-coupling approximation, in 
which the vortex dynamics is affected by the reservoir degrees of freedom only through 
the second order time correlation tensor (B(t)B), where the angular brackets indicate an 
average over the reservoir equilibrium ensemble and B(t) denotes a free time evolution 
for the reservoir operators. In addition, such tensor is naturally assumed to be isotropic 
in the x — y plane. We have also made use of the Markovian approximation which 
assumes that such correlations are short lived within observational times. From such 
master equation, we then extracted equations of motion for the expectation values of 
vortex position and momentum operators. Finally, after some algebra and by elimination 
of the momentum, we obtained the desired equation for (R(t)) f|. ||: 

m v (R) = i(m v n + y){R) . (6) 

The complex coefficient 7 is defined as 

7=?fV^T) e -^, (7) 



h Jo 

where ip(r) is the imaginary part of the time correlation isotropic tensor element, 

^(r)=Im(<5 j (r)5 j >) . (8) 

In the right hand side of equation (|6]) one can identify two forces, namely the 
Magnus force Q and the drag force. In particular, the drag force 

F D = i 1 (R) (9) 

has two components. One is parallel to the Magnus Force which we shall call the 
Transverse Force (TF) : 

Ft — ilZe(j) (R), (10) 

and the other one is parallel to the velocity and we shall refer to it as the Longitudinal 
Force (LF) : 

F L = -Im( 7 ) (R). (11) 

We want to recall that it is well known that the LF must be opposite to the velocity 
(in our notation this means Xm{pf) > 0). With respect to the TF great controversy still 
exists upon its direction |], ||, || or even more, upon its proper existence f7|, || . 

In the present work, the reservoir vector operator B of Eq. @ is chosen as to 
describe qp scattering by the vortex. In terms of creation a+ and annihilation a q 
operators for a qp of momentum q it reads, 

B = £(k-q)A fei(? a+a q , (12) 

k,q 
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being k — q the transferred momentum and Afc i9 a vortex-qp coupling constant which, in 
order to preserve isotropy, is assumed to depend only on the modulus of qp momentum 

(q = |q|)- 

To calculate the drag coefficients we must first compute the correlation function 
in (||). Taking into account the time evolution a q (£) = a q (0)e~™ 9 * and after some 
calculations, we have 

< Bj(t) Bj >= W\^ 2 (k 2 + q 2 )e^- w ^n(w k )[l + n(w q )} , (13) 

1 k,q 

where n(w q ) denotes the mean number of qp's of energy hw q , 

n(w q ) =< a+ a q >= _ - . (14) 



Note that at temperature T = the correlation function (|T^) vanishes due to the 
factor n(w k ), yielding a vanishing drag force as expected. This behavior would not be 
reproduced by a linear interaction in the qp operators, as for example the one proposed 
in the Hamiltonian of Ref . [f| . 

The parameter 7 (Eq. (0)) can be written in terms of the Fourier transform ip[w] 
of the imaginary part of the correlation function fllTf ) as 

2m( 7 ) = ^-^M (16) 
where V refers to the principal part and 

= 7 H( fc2 + q 2 )\^q\ 2 [n(w q ) - n(w k )]S(w k - w q - w) . (17) 

4 k,q 

The above expressions (|T6f) and (fL7|) for the dissipative coefficient in (|TTD can be 
interpreted from the scattering processes embodied in our interaction Hamiltonian. In 
fact, we first note that the vortex Hamiltonian ([I]) has an equally spaced level spectrum 
of separation TiVt (the so-called Landau levels |10|]) and, in addition, the operator v which 
couples the vortex to the reservoir qp's, can be expressed from (|j) as a linear combination 
of creation and destruction operators of a vortex energy quantum|lC|. Thus, according 
to @ and (|T2D, we may identify the scattering processes embodied in our model as those 
involving one vortex energy quantum HQ jointly with the creation and destruction of one 
qp. Such processes can also be identified from (|17]) as follows. Each term proportional 
to n(w q ) represents the process by which a qp of energy hw q combines with a vortex 
quantum HQ to form a qp of energy hw k = hw q + hQ (this may be seen from fliTf ) and the 
Dirac delta in (|TTD). The weight n{w q ) of this process can be easily understood by taking 
into account the thermal origin of the "incoming" qp of energy Tiw q in contrast to the 
interaction origin of the "outgoing" qp of energy hw k - Similarly, the terms weighted by 
n(wk) in ( |T7f ) represent the annihilation of a qp of energy hw k jointly with the creation 
of a vortex quantum hfl and a qp of energy hw q . Finally, we note that the positive sign 
of the dissipative coefficient Xm{pf) arises from the factor n{w q ) — n(w k ) > in ([H]) 
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i. e., the weight of processes involving vortex energy loss must be greater than those 
causing vortex energy gain. 

At low enough temperatures it is to be expected that only scattering processes 
involving phonons should be relevant. Therefore, we shall use the phonon dispersion 
relation w q = c s q (being c s the sound velocity) and impose a cutoff momentum to the 
qp's. This amounts to neglect all the scattering processes which yield not a phonon 
as the "outgoing" qp. Note that the "incoming" thermal qp will be surely a phonon. 
Under such an approximation Eq. (|17D can be written as, 



ip[w] = J dw' S(w', w' + w) \w' 2 + (w' + w) 2 } 

x [n(w') -n(w' + w)\ (18) 

where we introduce the so-called scattering function |IT| , defined as 

S(w', w") = £ \A k , q \ 2 5(w' - w q )5(w" - w k ) (19) 

k,q 

and being related to the scattering of the environmental excitations between states of 
frequencies w' and w". 

The integral in ( |TBD has to be solved numerically except in the limit T — > 0. In 
fact, for hw/ksT — > oo the gain term n{w' + w) can be neglected and the loss one n{w') 
"filters" all but the lowest frequency "incoming" phonons. This means that the factors 
accompanying n{w') in (|i~8l) can be approximated to lowest order in w' . In particular, 
the scattering function (fl9|) is assumed to be a continuous symmetric function of both 
variables [[□]] satisfying S(w',w) ~ S(w)w' p for w' — > 0. Thus equation ([18]) can be 



approximated for T — >• as follows, 



tp[w] ~ -w 2 S(w) I dw' w' p n(w') 







i , w . . s 2 of \ ( ^bT 



p+1 



= — p \((p+l) w iS(w)^-j , (20) 

where ((n) (n > 2) denotes the Riemann zeta function. 

It is convenient to notice that our model is unable to provide an a priori explicit 
form for the scattering function, because we treat vortex and qp excitations as separate 
entities which are assumed to interact by means of a generic Hamiltonian. Any 
additional information should be based upon experimental results or more fundamental 
theories. In fact, for the lowest temperature domain, only theoretical results are at 
present available and they predict a T 5 dependence for the LF [I]. Hence, we set p = 4 
in (|20| ) and accordingly a low frequency ~ w 4 behavior for the scattering function. To 
perform numerical calculations which illustrate our results we shall assume a simple form 



for the scattering function (|l^), as a generalization of the usual super-Ohmic dissipation 
with an exponential cutoff 



S(w', w") oc w' 4 e- w '/ w °w" 4 e- w " /w ° , (21) 
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being w Q a frequency cutoff parameter which allows us to select only the phonon part 
of the 4 He qp excitations spectrum i. e., the frequencies below w c ~1.2 ps" 1 . This 
is evident from Fig. 1 where we plot the one variable scattering function S(w,w) for 
w o =0.06 ps" 1 . In addition, we plot the frequency spectrum f(w) of the normal fluid 
density (i. e., p n = J °° f{w)dw) for T=0.4 K, which shows that at most the first half of 
the phonon spectrum makes a relevant contribution to the qp excitations according to 
a similar behavior for the one variable scattering function. 

From Eqs. (|15|) and (|l6j) we note that the value of cyclotron frequency, which in 
turn depends on the vortex mass m v , (see Eq. (f|)) is necessary for the determination 
of the drag force. Unfortunately, there is also a controversy regarding the calculation of 
the vortex mass|9|, [13|, [14], |15[ that makes the possible value of Q to range from a lower 
bound fill Q m in ~ 0.1 ps" 1 to the upper bound|§] Q ma x — 3 ps" 1 . This suggests that a 
study of the drag force dependence on Q could be illustrative. 

In Fig. 2 we show the coefficients of both, the transverse and the longitudinal 
components of the force, that is IZei^j) and Im('~f) respectively, as functions of Q for 
different temperatures. We see that, apart from a change of scale, the shape of curves 
does not change appreciably along the temperature range < T < OAK that we are 
considering. The only effect seems to be a shift to low frequencies with increasing 
temperatures (note that for a fixed vortex mass, the cyclotron frequency (|^) decrease 
with increasing temperature due to the factor p s , however this variation is negligible for 
T < OAK). 

Regarding Xm{pf) we first indicate that according to (|16|), the plots of this function 
in Fig. 2 turn out to be proportional to the Fourier transform of ip(t). In addition, we 
notice that Xm{pf) vanishes for Q > w c . Such absence of dissipation is simply understood 
as a direct consequence of energy conservation in the scattering processes. This is 
imposed by the Dirac delta in flT7[), since it reflects the fact that phonon-to-phonon 
scattering events can only take place if the vortex energy quantum does not exceed the 
most energetic phonon of the spectrum. Therefore, we may conclude that dissipation due 
to phonon^phonon scattering should be negligible unless the vortex mass yields a value 
of cyclotron frequency less than w c ~1.2 ps" 1 , which is an intermediate value between 
the quoted frequencies VL min and VL max . Moreover, the value %Vt maX) arising from a usual 
hydrodynamical prescription for the calculation of the vortex massfU, exceeds the energy 
of any undamped qp excitationfl^]. So, we may extend our previous conclusion by saying 
that any dissipation via qp-^qp scattering events should be negligible for such value of 
vortex mass, suggesting thus a possible scenario of dissipation via multi-qp scattering 
processes |L7] ]. 

Regarding the TF, we remark that if it were possible to draw some experimental 
information about at least the direction of this component, this would also shed more 
light on the value of the vortex mass. In fact, from Fig. 2 we note that at intermediate 
values between f2 min and w c , the TF changes its direction. In particular, a positive 
(negative) value of this component implies Q > 0.45 ps" 1 (f2 < 0.3 ps" 1 ), while a 
change of sign yields some intermediate value of Q. Unfortunately no experimental data 
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are available for temperatures below 0.4 K. 

It is important to note, in concluding this report, that our study has been restricted 
to a strictly rectilinear vortex, neglecting thus any possible contribution to the forces 
arising from vortex line curvatures. Clasically, vortex lines can be deformed as helical 
waves known as Kelvin waves, and there is experimental evidence of similar modes for 
quantized vortices in Helium II M. From the theoretical viewpoint, such waves could 
be regarded as a particular form of elementary excitation bound to the vortex. In fact, 
it has been shown that the elementary excitation spectrum of an imperfect Bose gas 
in presence of a straight vortex line consists of both, phonons and Kelvin-like waves 
[ i"8fl . Viewed from this perspective, being two different kinds of independent elementary 
excitations, the consequence of scattering processes involving phonons and Kelvin modes 
appears eventually as a possible secondary effect upon the values of longitudinal and 
transverse forces. This issue has been scarcely addressed in the literature. In particular, 
Sonin has long ago performed a simplified study [19j neglecting the vortex velocity due 
to Kelvin modes, but taking into account the effect of such modes on the scattering 
of phonons by means of an average over a classical Rayleigh- Jeans distribution for the 
oscillations of the vortex filament. He concluded that such effect could at most modify 
the calculated forces for a strictly rectilinear vortex in 1.7 % at T = 0.5K. Unfortunately, 
within our formalism a quantitative study of the dynamics of a true three-dimensional 
curved vortex filament, exhibits a high degree of difficulty 

In summary, starting from a microscopic model we have performed a calculation of 
the drag force on a moving vortex due to the scattering of qp excitations at temperatures 
below 0.4 K. We have also discussed, by analyzing both the longitudinal and transverse 
forces as functions of the cyclotron frequency, some important implications in the 
determination of possible scattering processes leading to dissipation, according to the 
values of vortex mass found in the literature. 
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Figure 1. One variable scattering function, S(w,w) for w o =0.06 ps -1 (dotted line) 
and frequency spectrum f(w) of the normal fluid density for T=0.4 K (solid line). 
Different scales were used for each function in order to normalize the size of both 
peaks. 
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Figure 2. Coefficients of the TF (TZe(jj) and the LF (Im( 7 )) for T ^ (Eq. 
(20)) (solid line); T = OAK (dashed line) and T = OAK (dotted line). Different 
scales were used for each temperature in order to normalize the size of the Tm{ r fj 
peak. The relative scale factors are 105.82 and infinity for T = OAK /T = OAK and 
T = OAK/T — > 0, respectively. The exponential parameter w o =0.06 ps _1 was used in 
Eq. (21). 



